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Abstract 

We report on double-differential inclusive cross-sections of the production of secondary 
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1 Introduction 

The HARP experiment arose from the realization that the inclusive differential cross-sections of 
hadron production in the interactions of low-momentum protons with nuclei were known only 
within a factor of two to three, while more precise cross-sections are in demand for several rea- 
sons. Consequently, the HARP detector was designed to carry out a programme of systematic 
and precise measurements of hadron production by protons and pions with momenta from 1.5 
to 15 GeV/c. It is shown schematically in Fig.[T] 




Fig. 1 : Schematic view of the HARP detector. 



The detector extended longitudinally over 14.7 m and combined a forward spectrometer 
with a large-angle spectrometer. The latter comprised a cylindrical Time Projection Chamber 
(TPC) around the target and an array of Resistive Plate Chambers (RPCs) that surrounded the 
TPC. The purpose of the TPC was track reconstruction and particle identification by dE/dx. 
The purpose of the RPCs was to complement the particle identification by time of flight. 

The HARP experiment was performed at the CERN Proton Synchrotron in 2001 and 2002 
with a set of stationary targets ranging from hydrogen to lead, including beryllium. 

We report on the large-angle production (polar angle 9 in the range 20° < 9 < 125°) of 
secondary protons and charged pions, and of deuterons and charged kaons, in the interactions 
with a 5% Aabs Be target of +8.9 GeV/c protons and pions, and of —8.0 GeV/c pions. 

The data analysis presented in this paper rests exclusively on the calibrations of the TPC 
and the RPCs that we, the HARP-CDP group, published in Refs. [1] and [2]. As discussed in 
Refs. [3] and [4], and succinctly summarized in this paper's Appendix, our calibrations disagree 
with calibrations published by the 'HARP Collaboration' [5-8]. Conclusions of independent 
review bodies on the discrepancies between our results and those from the HARP Collaboration 
can be found in Refs. [9, 10]. 
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2 The T9 proton and pion beams 

The protons and pions were delivered by the T9 beam line in the East Hall of CERN's Pro- 
ton Synchrotron. This beam line supports beam momenta between 1.5 and 15 GeV/c, with a 
momentum bite Ap/p ~ 1%. 

Beam particle identification was provided for by two threshold Cherenkov counters filled 
with nitrogen, and by time of flight over a flight path of 24.3 m. In the +8.9 GeV/c and 
—8.0 GeV/c beams, the pressure of the nitrogen gas was set such that protons were below 
threshold for Cherenkov light but pions above. The time of flight of each beam particle was 
measured by three scintillation counters with a precision of 106 ps^^. 

Figure |2] (a) shows the relative velocity (3 from the beam time of flight of positive parti- 
cles in the +8.9 GeV/c beam, with protons distinguished from 'pions'^^ by the absence of a 
beam Cherenkov signal. Vice versa, Fig.[2](b) shows the signal charge in one beam Cherenkov 
counter, with protons and pions distinguished by the signal charge in the other beam Cherenkov 
counter. All measurements are independent of each other and together permit a clean separation 
between protons and pions, respectively, with a negligible contamination of less than 0.1% by 
the other particle. 
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Fig. 2: (a) Relative velocity [3 from the beam time-of-flight system of protons (shaded histogram) and 
pions, with the particles identified in the beam Cherenkov counters; (b) charge response of the light 
signal in one beam Cherenkov counter from protons (shaded histogram) and pions, with the particles 
identified in the other beam Cherenkov counter. 



The pion beam had a contamination by muons from pion decays. This contamination was 
measured to be (1.7 ± 0.5)% of the pion component of the +8.9 GeV/c beam [11]. For the 
—8.0 GeV/c beam, this contamination is (1.9 ± 0.5)%. The pion beam also had a contamination 
by electrons from converted photons from 7r° decays. This contamination was determined to be 

^^Under stable conditions of the beam optics, such that an average particle velocity could be used, the time-of- 
flight precision could be improved to 77 ps; in our analysis, no use was made of this option. 

^^The 'pions' comprise small contaminations by muons and electrons, indistinguishable both by time of flight 
and by beam Cherenkov signals. 
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(1.2 ± 0.5)% of the pion component of the +8.9 GeV/c beam [12]. We take the same electron 
fraction for the —8.0 GeV/c beam. For the determination of interaction cross-sections of pions, 
the muon and electron contaminations must be subtracted from the incoming flux of pion-like 
particles. 

The beam trajectory was determined by a set of three multiwire proportional chambers 
(MWPCs), located upstream of the target, several metres apart. The transverse error of the 
projected impact point on the target was 0.5 mm from the resolution of the MWPCs, plus a 
contribution from multiple scattering of the beam particles in various materials. Excluding the 
target itself, the latter contribution is 0.2 mm for a 8.9 GeV/c beam particle. 

The size of the beam spot at the position of the target was several millimetres in diameter, 
determined by the setting of the beam optics and by multiple scattering. The nominal beam 
position^) was at Xbcam = 1/beam = 0, howcvcr, excursions by several millimetres could occur^^ 
A loose fiducial cut ^J a^beam + 1/boam < mm ensured full beam acceptance. The muon and 
electron contaminations of the pion beam, stated above, refer to this acceptance cut. 

We select 'good' beam particles by requiring the unambiguous reconstruction of the particle 
trajectory with good x^- In addition we require that the particle type is unambiguously identi- 
fied. We select 'good' accelerator spills by requiring minimal intensity and a 'smooth' variation 
of beam intensity across the 400 ms long spill^\ 

3 The large-angle spectrometer 

In harp's large-angle region, a cylindrical TPC [1] had been chosen as tracking detector. It 
was embedded in a solenoidal magnet that generated a magnetic field of 0.7 T parallel to the 
TPC axis. The magnet was in general operated with its polarity tied to the beam polarity 

The TPC filled most of the inner bore of the magnet, leaving a 25 mm wide gap between 
TPC and magnet coils. This gap was used to house two overlapping layers of 2 m long RPCs [2] 
directly mounted onto the outer field cage of the TPC. 

The layout of the TPC and its position in the solenoidal magnet is shown in Fig. [3} The TPC 
has an external diameter of 832 mm and an overall length of ~2 m. It consists of two Stesaht 
cylinders forming the inner and outer field cages, a wire chamber with pad readout, located at 
the upstream end, and a high- voltage (HV) membrane at 1567 mm distance from the pad plane. 
The inner field cage extends over about half of the drift volume; it encloses the target, the centre 
of which is located 500 mm downstream of the pad plane. 

The tracking volume extends radially from 75 mm to 385 mm and over ~1.5 m longitu- 
dinally. Electrons from ionization induced by charged particles in the TPC gas drift upstream 
under the influence of the longitudinal electrical field; they are amplified in the wire chamber 
and read out through pads arranged in six identical sectors, as shown in Fig. |3| Each sector 
comprised 662 readout pads of dimensions 6.5 x 15 mm^ arranged in 20 concentrical rows. 

Our calibration work on the HARP TPC and RPCs is described in Refs. [1] and [2], and in 
references cited therein. In particular, we recall that static and dynamic TPC track distortions 

■^^A right-handed Cartesian and/or spherical polar coordinate system is employed; the z axis coincides with the 
beam line, with +z pointing downstream; the coordinate origin is at the centre of the beryllium target, 500 mm 
downstream of the TPC's pad plane; looking downstream, the +a: coordinate points to the left and the +j/ coordi- 
nate points up; the polar angle Q is the angle with respect to the +z axis; when looking downstream, the azimuthal 
angle </> increases in the clockwise direction, with the +2; axis at = 0. 

■^^The only relevant issue is that the trajectory of each individual beam particle is known, whether shifted or not, 
and therefore the amount of matter to be traversed by the secondary hadrons. 

^^A smooth variation of beam intensity eases corrections for dynamic TPC track distortions. 

^^The HARP data-taking convention was that i?z > refers to positive beam polarity. 
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Fig. 3: Longitudinal cut through the TPC and the solenoidal magnet; the beam enters from the left side; 
the small figure to the right shows the layout of the six TPC readout sectors, looking downstream. 



up to ~10 mm have been corrected to better than 300 /xm. Therefore, TPC track distortions do 
not affect the precision of our cross-section measurements. 

4 Track reconstruction and particle identification 

4.1 Pattern recognition in the TPC 

The clusters measured by the TPC constitute space points along the track trajectory. Each space 
point has three uniquely determined coordinates: r, 0, and z. Our pattern recognition of tracks 
with > 0.05 GeV/c originating from the target region is based on the TOPAZ histogram 
technique [13]: a 2-dimensional histogram of the ratio zjr against azimuthal angle is filled 
with all reconstructed clusters. Physical tracks populate one or two adjacent bins (the bin sizes 
are suitably chosen) and thus are easily recognised. 

4.2 Helix fit of TPC tracks 

For the fit of trajectories in the TPC we adopted the 'Generalized Least-Squares Fit' (GLSF) 
concept. This is the formal generalization of the standard least-squares fit for an arbitrary num- 
ber of error-prone dimensions, and the solution of the equations resulting from the -y^ mini- 
mization with the Lagrange-multiplier method. The mathematical intricacies can be found in 
Ref. [14]. For the three parameters that describe the circle projection of a helix, we adopted 
the TOPAZ parametrization [15], for the attractive feature of avoiding any discontinuity in the 
numerical values of fit parameters. Most importantly, it features a smooth transition between 
charge signs of a track. For more details on the parametrization and the fit procedure, we refer 
to Ref. [16]. 

The GLSF must start from reasonable starting values of the parameters that describe the 
helix. They are obtained by the Chernov-Ososkov least-squares algorithm [17]. 

Our GLSF algorithm yields the transverse momentum of a track, its charge sign, its polar 
angle 6*, and its closest point of approach to the z axis. 
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4.3 Virtual beam point 

The pt resolution of tracks can be significantly improved by the use of the beam point''^ as 
an additional point to the trajectory in the TPC. The transverse coordinates of the beam point 
are known from the extrapolation of the trajectory of the incoming beam particle. Their errors 
originate from three sources. The first is from the extrapolation error of the beam trajectory that 
is measured by MWPCs; the second stems from multiple scattering of the beam particle; and 
the third from multiple scattering of the secondary particle in materials between the vertex and 
the TPC volume. 

However, the correct error assignment to the beam point is not sufficient. Since a secondary 
track loses energy by ionization in the target and in materials between the vertex and the TPC 
volume, a correction must be calculated that replaces the real beam point by a 'virtual' beam 
point which is bias-free with respect to the extrapolation of the trajectory measured in the TPC. 
It is this virtual beam point, and not the real beam point, that is used in the (final) track fit. It is 
determined in an iterative procedure that starts from the fit of the track momentum in the TPC 
gas, including the real beam point. The fitted trajectory in the TPC gas is then back-tracked 
to the beam particle trajectory taking the energy loss and multiple scattering into account. It 
renders a first estimate of the virtual beam point. Using this estimate, the track in the TPC 
gas is again fitted, and the procedure is iterated until the position of the virtual beam point 
is stable. Since in the calculation of the move from the real to the virtual beam point the 
energy loss is taken into account, and since the energy loss depends on the type of particle, 
three different virtual beam points are calculated according to the proton, pion, and electron 
hypotheses. Accordingly, three different track fits are performed. 

The fit with the virtual beam point included gives the best possible estimate of the particle 
momentum in the TPC gas. In order to determine what is really needed, namely the momentum 
at the vertex, in a last step the particle is tracked back to the vertex, taking into account the 
energy loss under the three different particle hypotheses. The track parameters at the vertex are 
used for the determination of differential cross-sections. 

4.4 Particle identification algorithm 

The particles detected in HARP's large-angle region are protons, charged pions, and electrons®-' 
(we disregard here small admixtures of kaons and deuterons which will be discussed in Sec- 
tion |7]). The charged pion sample comprises muons from pion decay since the available instru- 
mentation does not distinguish them from charged pions. 

The dE/dx and the time-of-flight methods of particle identification are considered indepen- 
dent. 

To separate measured particles into species, we assign to each particle a probability of being 
a proton, a pion (muon), or an electron, respectively. The probabilities add up to unity, so that 
the number of particles is conserved. 

Each track is characterized by four measured quantities: (transverse momentum), 9 (po- 
lar angle), P (relative velocity) and dE/dx (specific ionization). For particle identification 
purposes, these variables refer to reconstructed ('smeared') variables in both the data and the 
Monte Carlo simulation. 

In every bin of (p^, 9), the probability V{i\P, dE/dx, pt, 9) of a particle to belong to species 
i{i = \ [proton], 2 [pion], 3 [electron]) in a mixture of protons, pions, and electrons is according 

^^The 'beam point' is the best estimate of the interaction vertex of the incoming beam particle. 
*^The term 'electron' also refers to positrons. 
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to Bayes' theorem as follows: 



V{i\(3,dE/dx,pT:,e) = — — — , (1) 

X;i=i [P\I3, dE/dx|z, pt, 0) ■ P{h Pt, d)] 

where the sum 

3 

j2n^\p^dE/dx,PT,e) 

i=l 

is normalized to unity. The probabilities P{i,pT, 0) are given by 



P{i,PT, 



where Ni{pT,9) is the number of particles of species i in the respective data sample. Then 
Eq. ([T]) becomes 

-pri/^^p/^ m PiP,dE/dx\t,pT,e)-N,{p^,e) 
Vn\p,dE/dx,pT,0) = — ^ . (2) 

j:-[mdE/dx\t,pT,e)-N,{p^,9)] 

We note that in Eqs. ([T]) and ^ the term P{P, dE /dx\i,pT, 9) denotes a probability density 
function which is normalized to unity. This probability density function must represent the data 
in the bin (px, d). 

Before determining the probability represented by Eq. Q, the probability density functions 
P{P,dE/dx\i,pT,9) and the particle abundances Ni(pT,9) must be known. This seemingly 
circular situation is resolved by an iterative comparison of data with the Monte Carlo simulation, 
to achieve agreement of the distributions in both variables P and dE/dx. With a view to starting 
from abundances as realistic as possible, the comparison is initially limited to regions in phase 
space where the particle species are unambiguously separated from each other in either dE / dx 
or p. In other words, the few parameters that govern the probability density functions and the 
particle abundances are determined from the data in every bin of (pr, 9). 

In case one of the two identification variables is absent^^, only the other is used. In the rare 
cases where both identification variables are absent, the identification probabilities reproduce 
the estimated particle abundances. 

4.5 Particle abundances 

Particle abundances cannot a priori be expected to be correct in the Monte Carlo simulation. 
Therefore in general the particles must be weighted such that data and Monte Carlo distributions 
agree. 

We had expected that the Monte Carlo simulation tool kit Geant4 [18] would provide us 
with reasonably realistic spectra of secondary hadrons. We found this expectation more or less 
met by Geant4's so-called QGSP_BIC physics list, but only for the secondaries from incoming 
beam protons. For the secondaries from incoming beam pions, we found the standard physics 
lists of Geant4 unsuitable [19]. 

To overcome this problem, we built our own HARP CD? physics list for the production of 
secondaries from incoming beam pions. It starts from Geant4's standard QBBC physics list, 
but the Quark-Gluon String Model is replaced by the FRITIOF string fragmentation model for 

^^For example, because of too few clusters to calculate clE/dx, or a missing RPC pad. 
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kinetic energy E > Q GeV; for < 6 GeV, the Bertini Cascade is used for pions, and the 
Binary Cascade for protons; elastic and quasi-elastic scattering is disabled. Examples of the 
good performance of the HARP_CDP physics list are given in Ref. [19]. 

Figure |4] demonstrates the level of overall agreement between data and Monte Carlo simula- 
tion in the variable l/pr, after convergence of the iterative procedure to determine the smooth 
weighting functions to the latter. The figure also shows, for incoming protons and for a typi- 
cal polar-angle range, the subdivision of the data into particle species by applying the particle 
identification weights. 




Fig. 4: 1/pt spectra of the secondary particles from +8.9 GeV/c beam protons on a 5% Aabs Be target, 
for polar angles 50° < 9 < 60°; black triangles denote data, the solid lines Monte Carlo simulation; 
the shaded histograms show the subdivision of the data into particle species by applying the particle 
identification weights: light shading denotes protons, medium shading pions, and dark shading electrons. 

Once the abundances are determined, for any pair of dE/dx and P, and using the experi- 
mental resolution functions, the probability can be derived that the particle is a proton, a pion, 
or an electron. This probability is consistently used for weighting when entering tracks into 
plots or tables. 

5 Physics performance 

5.1 Physics performance of the TPC 

From the requirement that a tt"*" and a vr" with the same RFC time of flight have the same 
momentum, and from the error of the magnetic field strength which is less than 1%, the absolute 
momentum scale is determined to be correct to better than 2%, both for positively and negatively 
charged particles. 
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Figure |5] (a) shows the l/px difference for positive particles with 0.6 < /? < 0.75 and 
45° < 9 < 65°, between the measurement in the TPC and the determination from RPC time 
of flight with the proton-mass hypothesis. The selection cuts ensure a practically pure sample 
of protons (the background from pions and kaons is negligible as suggested by the very small 
contribution of negative particles selected with the same cuts that are shown as dots in Fig.[5](a)). 
A net TPC resolution of a^l/pi) = 0.20 (GeV/c)^^ is obtained by subtracting the contribution 
of ~0.18 (GeV/c)^^ from the time-of-flight resolution and fluctuations from energy loss and 
multiple scattering in materials between the vertex and the TPC volume quadratically from the 
convoluted resolution of 0.27 (GeV/c)~^. 

Figure|5](b) shows the net TPC resolution o"(1/pt) as a function of p. Figure[5](c) shows the 
same as a function of 9. The agreement with the expectation from a Monte Carlo simulation is 
satisfactory. The resolution a(l/p-p) is typically 20% and worsens towards small P and small 9. 
This is because in both cases the position error of the virtual beam point increases owing to 
increased multiple scattering in materials before the protons enter the TPC. 




Fig. 5: (a) Difference of the inverse transverse momenta of positive (shaded histogram) and negative 
(dots) particles from the measurement in the TPC and from the determination from RPC time of flight, 
for 0.6 < j3 < 0.75 and for 45° < 6 < 65°; the positive particles are protons, with a very small 
background from pions and kaons; (b) (j(1/pt) of protons with 45° < 6 < 65° as a function of their 
relative velocity /3; (c) c7(1/pt) of protons with 0.6 < /3 < 0.75 as a function of their polar angle 6. 



Data from the elastic scattering of incoming pions or protons on protons at rest have the 
added feature of a kinematical constraint. The possibility to calculate from the four-momentum 
of the incoming beam particle and the polar angle 9 the momentum of the large-angle recoil 
proton, permits a valuable cross-check of the TPC's pt resolution. Figure |6] shows the result 
from the elastic scattering of incoming +3 GeV/c protons and vr+'s in a liquid hydrogen target. 
Here, the of the recoil proton has been determined in the following two ways: l/p'^^^ is 
determined from the reconstructed track curvature in the TPC; 1/p^^'^ is predicted from the 
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elastic scattering kinematics from the polar angle of the recoil proton which is little affected 
by TPC track distortions. Figure |6] demonstrates a resolution of cr{l/p^) ~ 0.19 (GeV/c)^^ 
after unfolding a contribution of a{l /pt) ~ 0.12 (GeV/c)~^ from fluctuations from energy loss 
and multiple scattering in materials between the vertex and the TPC volume. The measured 
difference in is 0.8%, in line with the 2% uncertainty of the momentum scale. 




1/p^(nneas) - 1/p^(pred) [GeV/c] 



Fig. 6: The difference 1 /p^'^^^ — 1 /p^'^'^ from large-angle recoil protons in elastic scattering events from 
a +3 GeV/c beam impinging on a liquid hydrogen target; the tail at the right side reflects the Landau tail 
in the proton energy loss in materials between the interaction point and the TPC volume. 



The polar angle 6 is measured in the TPC with a resolution of ~9 mrad, for a representative 
angle of 6* = 60°. To this a multiple scattering error has to be added which is ~7 mrad for a 
proton with = 500 MeV/c and 9 = 60°, and ~4 mrad for a pion with the same characteristics. 
The polar-angle scale is correct to better than 2 mrad. 

Besides the and the polar angle 9 of tracks, the TPC also measures dE/dx with a view 
to particle identification. The dE/dx resolution is 16% for a track length of 300 mm. 

5.2 Physics performance of the RPCs 

The intrinsic efficiency of the RPCs that surround the TPC is better than 98%. While the system 
efficiency for pions with p^ > 100 MeV/c at the vertex is close to the intrinsic efficiency, it is 
slightly worse for protons because of their higher energy loss in structural materials. Protons 
with p < 350 MeV/c at the vertex get absorbed before they reach the RPCs and thus escape 
time-of-flight measurement. 

The intrinsic time resolution of the RPCs is 127 ps and the system time-of-flight resolution 
(that includes the jitter of the arrival time of the beam particle at the target) is 175 ps. 

Figure [t] (a) shows the specific ionization dE/dx, measured by the TPC, and Fig. |7](b) 
the relative velocity P from the RPC time of flight, of positive and negative secondaries, as 
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a function of the momentum measured in the TPC. The figures demonstrate that in general 
protons and pions are well separated. They also underline the importance of the complementary 
separation by RPC time of flight at large particle momentum. 
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Fig. 7: (a) Specific ionization dE/dx and (b) velocity /3, versus the charge-signed momentum of positive 
and negative tracks in +8.9 GeV/c data; the boxes at the right side indicate the event statistics. 



6 Normalized secondary particle flux 

The measurement of the inclusive double-differential cross-section d^cr/dpdf2 requires the flux 
of incoming beam particles, the number of target nuclei, and the number of secondary particles 
in bins of momentum p and polar angle 9. We shall discuss these elements in turn. 



6.1 Beam intensity 

The event trigger had two levels. A first, loose, level required only time-coincident signals from 
beam scintillation counters. Irrespective of an interaction in the target, each 64th coincidence 
signal requested data readout as 'beam trigger' . A second, tighter level required in addition a 
signal in a cylindrical scintillator fibre detector that surrounded the target region, or a signal in a 
plane of scintillators in the forward direction (termed 'FTP' in Fig.|3]). Each such 'event trigger' 
also requested data readout. 

To achieve the wanted event statistics, the experiment was typically run with a dead time 
in excess of 50%, given the 400 ms long accelerator spill and a readout time of order 1 ms 
per event. Since the dead time affects the beam trigger and the event trigger in the same way, 
it cancels in the cross-section calculation. For a given data set, the flux of incoming beam 
particles is defined by the number of beam triggers, multiplied by the scale-down factor of 64. 
It is imperative, though, that the same cuts on the quality of the trajectory of the beam particle 
and on its identification be applied for accepted beam triggers and for accepted event triggers. 

The efficiencies of both the beam trigger and the event trigger are very close to 100%, 
thanks to majority requirements. The beam trigger efficiency cancels. For the event trigger, we 
determined an efficiency of (99.0±0.2)%. 
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6.2 Target 

The target was a cylinder made of high-purity (99.95%) beryllium, with a density of 1 .85 g/cm^, 
a radius of 15 mm, and a thickness of 20.5 ± 0.1 mm (5% Aabs)- 

The finite thickness of the target leads to a small attenuation of the number of incident beam 
particles. The attenuation factor is /att = 0.975. 

6.3 Track counting in bins of p-r and 6 

This paper is concerned with determining inclusive cross-sections of secondaries from the in- 
teractions of protons and pions with beryllium nuclei. This means that for a given data set, the 
secondaries are weighted with their probability of being a proton, a pion, or an electron, counted 
in bins of and 9, and related to the number of incoming beam particles and the number of tar- 
get nuclei. The counting of secondaries is done in an integral way without regard to track-event 
relations. 

Electrons stem primarily from the conversion of photons from vr" decays. They tend to 
concentrate at small momenta. Below 150 MeV/c, they are identified by both dE/dx and 
time of flight from the RPCs. From 150 to 250 MeV/c, the dE/dx of pions and electrons 
coincides and they are only identified by time of flight. The Geant4 electron abundance is 
compared with data in the region of good separation, as a function of momentum, and weighted 
to agree with the data. In the region of bad separation the electrons are subtracted using not 
the electron abundance predicted by Geant4, but the weighted prediction extrapolated from the 
region of good separation. Therefore, the Geant4 prediction is used only through its extrapolated 
prediction of the energy dependence of electrons with momentum larger than 250 MeV/c. 

Since the particle identification algorithm assigns to every particle a probability of being 
a proton, a pion, or an electron, the elimination of electrons from the samples of secondary 
protons and pions is straightforward. 

It is justified to think of secondary tracks as originating exclusively from proton and pion 
interactions: interactions of beam electrons might occasionally lead to low-momentum electron 
or positron tracks in the TPC, however, such tracks are recognized by the particle identification 
algorithm and disregarded in hadron production cross-sections. Interactions of beam muons can 
be neglected. 

Kaon and deuteron secondaries are initially part of pions and protons, respectively. Their 
identification is dealt with in Section |7l 

6.4 Track selection cuts 

We have a selection of 'good' TPC sectors: we discard tracks from the 'horizontal' sectors 2 
and 5 out of the six sectors (see Fig. |3]) for reasons of much worse than average performance, 
and the lack of reliable track distortion corrections [1]. 

Tracks are accepted if there are at least 10 TPC clusters along the trajectory. 

A cut in the azimuthal angle (f) is applied to avoid the dead regions of the six 'spokes' that 
subdivide the TPC pad plane into six sectors: 10° on one side and 2° on the other side of each 
spoke for tracks of one charge, and vice versa for the other charge. The asymmetric cut is 
motivated by the opposite bending of positive and negative tracks in the magnetic field. 

The polar- angle range of tracks is limited to the range 20° < 9 < 125°. Tracks are also 
required to point back to the target, within the resolution limits. 
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6.5 Correction for inefficiencies of track reconstruction and track selection 

The track reconstruction efficiency was determined by eyeball scanning of several thousand 
events by several physicists, with consistent results among them. The large number of scanned 
events permits us to determine the reconstruction efficiency as a function of geometric or kine- 
matic variables. For example, Fig. [8] (a) shows the reconstruction efficiency as a function of 
1/j9t for all cases where the human eye finds at least five (out of a maximum of 20) clusters 
along a trajectory. The average reconstruction efficiency is between 95% and 91%^^\ where the 
2% range reflects the variation between different data sets. 

We cross-checked the track reconstruction efficiency by requiring an RPC hit and at least 
two TPC clusters in the cone that is subtended by the respective RPC pad, as seen from the 
vertex. Figure [8] (b) shows the resulting reconstruction efficiency as a function of the track's 
azimuthal angle. Outside the TPC spokes and within the four 'good' TPC sectors, the recon- 
struction efficiency determined that way agrees with the result from the eyeball scan. 

The requirement of a minimum of 10 TPC clusters per track entails a loss that must be 
accounted for. Since the TPC cluster charge is in general larger for protons^^^ than for pions, 
the loss from this cut is different for protons and pions. Figure [8] (c) shows the efficiency of 
requiring 10 or more TPC clusters as a function of 1/pt, separately for protons and pions (the 
average number of clusters was ~14). 

The overall track efficiency was taken as the product of the track reconstruction efficiency 
and the probability of having at least 10 clusters along the trajectory. 
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Fig. 8: (a) Track reconstruction efficiency from the eyeball scan; (b) track reconstruction efficiency from 
RPC hits; the vertical lines denote TPC sector boundaries where the TPC 'spokes' render the efficiency 
lower; the azimuthal ranges of the two TPC sectors not used for the analysis are indicated by arrows; (c) 
efficiency of 10 or more TPC clusters for protons (open circles) and for pions (black circles); the errors 
shown are statistical only and mostly smaller than the symbol size. 



6.6 Further corrections 

In this section, we discuss a few more corrections that are applied to the data. In general, they 
are determined from a Monte Carlo simulation that reproduces the migration of track param- 
eters from generated ('true') to reconstructed ('smeared') ones. This concerns effects arising 

^^^This average holds for tracks with \pt\ > 0.1 GeV/c and 20° < < 125°, with safe distance from the 
insensitive azimuthal regions caused by the TPC 'spokes'. 
^^^The lower proton velocity leads to higher specific ionization. 
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from finite resolution, charge misidentification, pion decays into muons, and re-interactions of 
secondaries in materials between the vertex and the TPC volume^^^. There is also backscatter- 
ing of particles from the solenoid coil at large radius back into the TPC, however, tracks from 
backscattering are eliminated by the requirement that they originate from the target. 

Other than for the transverse momentum p^, migration is nearly negligible in the measure- 
ment of the polar angle 9. 

Charge misidentification occurs only at large transverse momentum, at the level of a few per 
cent. For example, a few 'antiprotons' at large transverse momentum are charge-misidentified 
protons and treated accordingly in the migration correction. 

Pion decay into muons occurs at the typical level of 2%. When the pion decay occurs in the 
first few centimetres of the flight path, the phenomenon is taken care of by the migration correc- 
tion. When the pion decay occurs later, the track is likely to be lost because of the requirement 
that it originates from the target. Therefore, each pion receives a weight that compensates on 
the average for the loss from decay along a path of 200 mm length. 

The re-interaction of secondaries takes place in the target material or in other materials 
between the target and the TPC volume. The typical probability for re-interaction is 3% for the 
former, and 2% for the latter. The re-interaction leads to tracks with other parameters than the 
initial track, and is taken care of by the migration correction. 

6.7 Systematic errors 

The systematic precision of our inclusive cross-sections is at the few-per-cent level, from er- 
rors in the normalization, in the momentum measurement, in particle identification, and in the 
corrections applied to the data. 

The systematic error of the absolute flux normalization is taken as 2%. This error arises 
from uncertainties in the target thickness, in the contribution of large-angle scattering of beam 
particles, in the attenuation of beam particles in the target, and in the subtraction of the muon 
and electron contaminations. Another contribution comes from the removal of events with an 
abnormally large number of TPC hits above threshold. 

The systematic error of the track finding efficiency is taken as 1 % which reflects differences 
between results from different persons who conducted eyeball scans. We also take the statistical 
errors of the parameters of a fit to scan results as shown in Fig. [8] (a) as systematic error into 
account. The systematic error of the correction for losses from the requirement of at least 10 
TPC clusters per track is taken as 20% of the correction which itself is in the range of 5 to 
30%. This estimate arose from differences between the four TPC sectors that were used in our 
analysis, and from the observed variations with time. 

The systematic error of the scale is taken as 2% as discussed in Ref. [1]. 

The systematic errors of the proton, pion, and electron abundances are taken as 10%. We 
stress that errors on abundances only lead to cross-section errors in case of a strong overlap of 
the resolution functions of both identification variables, dE/dx and p. The systematic error of 
the correction for migration, absorption of secondary protons and pions in materials, and for 
pion decay into muons, is taken as 20% of the correction, or 1% of the cross-section, whichever 
is larger. These estimates reflect our experience with remanent differences between data and 
Monte Carlo simulations after weighting Monte Carlo events with smooth functions with a 
view to reproducing the data simultaneously in several variables in the best possible way. 

^^'These re-interactions, especially in the target material, are different from re-interactions of secondaries in the 
nuclear matter of the same nucleus in which the incoming beam particle interacted; the latter is an integral part of 
the inclusive cross-section reported in this paper. 
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All systematic errors are propagated into the momentum spectra of secondaries and then 
added in quadrature. 

7 Kaon and deuteron production 

The statistics from the +8.9 GeV/c beam on a 5% Aabs beryllium target is much larger than for 
any other combination of beam and target. This permits us to investigate in this particular data 
set the production of K+'s and deuterons in addition to the dominant protons, vr+'s, and vr^'s. 
With a view to benefiting from the cancellation of systematic errors, we present results in terms 
of the ratios K+/7r+ and d/p. 

7.1 Kaons 

Figure |9] shows the relative velocity (3 of positive secondaries for the polar-angle range 20.5° < 
6 < 25.3° and momentum between 520 and 560 MeV/c. A logarithmic scale is employed to 
make K+ production visible which is at the level of a few per cent of the 7r+ production. The 
K"*" signal shows up between the proton and vr^ signal thanks to the good resolution of the (3 
measurement by the RPCs. The K+ signal is fitted with a Gaussian. The 7r+ signal is represented 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 

Particle |3 = v/c 



Fig. 9: Distribution of (3 of positive secondaries on a logarithmic scale, with the K+ signal showing up 
between pions and protons; the crosses show the f3 distribution of vr^'s; the shown fit is explained in the 
text. 



by a Gaussian together with a tail that is experimentally determined from the (3 distribution of 
the TT^'s. The latter is shown with crosses in Fig.|9j A possible contribution is minimized 
by a dE/dx cut. 

In order to maximize the time of flight and hence the separation power, we restrict the 
analysis to the forward region in the range 20° < 9 < 32°. The momentum is required to be in 
the range 400 < p < 700 Me V/c, and dE/dx must be between 70% and 155% of the nominal 
value. 
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Several corrections must be made to the fit results of the relative K+ abundance. Correcting 
for cuts on the charge of the RPC signal, made with a view to optimizing time-of-flight reso- 
lution, reduces the signal by 5%. The correction for the non-Gaussian tail of the (3 distribution 
of K+'s increases the signal by 8%. The correction for different absorption of K+'s and vr+'s in 
structural materials increases the signal by 1%. 

Altogether, the resulting ratio is 

averaged over the said range of momentum and polar angle, and over the proton and 7r+ beams 
(in kaon production, no significant difference is seen between these beams). Figure [10] shows 
the K+/7r+ ratio as a function of particle momentum and compares the measured ratios with 
the ratios from the FRITIOF and Binary Cascade hadron production models in Geant4. The 
data points are closer to the prediction by the FRITIOF model, however, the dependence on 
momentum is not reproduced. The agreement with the Binary Cascade model is poor. 




0.4 0.5 0.6 0.7 0.8 

Momentum at vertex (GeV/c) 



0.9 



Fig. 10: Ratio of K"'"/7r"'" in +8.9 GeV/c proton and 7r+ interactions with beryllium nuclei; the data points 
(black circles) are closer to the prediction by the FRITIOF model (open histogram) in Geant4, but agree 
poorly with its Binary Cascade model (shaded histogram). 



7.2 Deuterons 



Figure 1 1 shows the dE / dx of positive secondaries for the polar-angle range 30° < 9 < 45° and 
the momentum range from 500 to 600 MeV/c (this momentum range refers to the momentum 
measured in the TPC and not to the momentum at the vertex). Pions and electrons are reduced 
by a loose time-of-flight cut. A clear signal of deuterons is visible at large dE/dx, next to the 
abundant protons. 

In order to transform the ratio measured in the TPC volume to that at the vertex, appropriate 
corrections for the different energy loss of protons and deuterons in materials between the vertex 
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log(dE/dx) 

Fig. 1 1 : Distribution of dE / dx of positive secondaries, with the deuteron signal showing up at large 
d£'/dx; see the text for the cuts applied; both the proton and deuteron signals are fitted with a Gaussian. 



and the TPC volume, and for differences in the momentum spectra of protons and deuterons, 
must be applied. The results for the ratio 

averaged over the momentum at the vertex between 600 MeV/c and 1050 MeV/c, are given in 
Tabled 

Table 1: Ratio of deuterons to protons, for different beam particles, averaged over the momentum at 
the vertex between 600 MeV/c and 1050 MeV/c. 





+8.9 GeV/c protons 


+8.9 GeV/c 71+ 


-8.0 GeV/c n- 


20° < ^ < 30° 
30° < ^ < 45° 
A5° <e < 65° 
65° <9 < 90° 
90° <e < 125° 


0.051 ±0.004 
0.076 ± 0.004 
0.113 ±0.009 
0.17 ±0.03 
0.26 ±0.04 


0.038 ±0.004 
0.065 + 0.005 
0.094 ±0.008 
0.15 ±0.03 
0.21 ±0.05 


0.05 ±0.01 
0.07 ±0.01 
0.13 ±0.02 
0.16 ±0.03 
0.30 ±0.05 



The ratios for the +8.9 GeV/c proton and 7r+ beams are shown in Fig. 12 for the polar- 
angle range 30° < 9 < 45° as a function of the momentum at the vertex. We note that the 
deuteron abundance is resonably well reproduced by the FRITIOF String Fragmentation model 
used in the Geant4 simulation tool kit, while it is underestimated by about one order of magni- 
tude by the Binary Cascade model. 
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Fig. 12: Ratio in +8.9 GeV/c proton and 7r+ interactions with beryllium nuclei, for 30° < 6 < 45°, 
as a function of the momentum at the vertex; black data points refer to the proton beam, open circles to 
the 7r+ beam. The proton beam data are compared with the predictions of the FRITIOF (open histogram) 
and the Binary Cascade (shaded histogram) models in Geant4. 



8 Double-differential inclusive cross-sections of protons and pions 
In Tables |2}jl0 we give the double-differential inclusive cross-sections d'^a/dpdQ for all nine 
combinations of incoming beam particle and secondary particle, including statistical and sys- 
tematic errors. In each bin, the average momentum and the average polar angle are also given. 

Cross-sections are only given if the total error is not larger than the cross-section itself. 
Since our track reconstruction algorithm is optimized for tracks with above ~70 MeV/c in 
the TPC volume, we do not give cross-sections from tracks with below this value. Because 
of the absorption of slow protons in the material between the vertex and the TPC gas, and with 
a view to keeping the correction for absorption losses below 30%, cross-sections from protons 
are limited to p > 350 MeV/c at the interaction vertex. Proton cross-sections are also not given 
if a 10% error on the proton energy loss in materials between the interaction vertex and the TPC 
volume leads to a momentum change larger than 2%. Pion cross-sections are not given if pions 
are separated from protons by less than twice the time-of-flight resolution. 

The data given in Tables [2- 10 are available in ASCII format in Ref. [20]. 
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Table 2: Double-differential inclusive cross-section d^a/dpdQ [mb/(GeV/c sr)] of the production of 
protons in p -I- Be ^ p -I- X interactions with +8.9 GeV/c beam momentum; the first error is statistical, 
the second systematic; pt in GeV/c, polar angle 9 in degrees. 
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0.459 


34.9 


33.39 


± 
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± 
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± 
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± 
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0.42-0.50 
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0.18 
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0.11 


± 
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96.5 
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0.13 


± 


0.37 
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2.28 


± 


0.07 


± 


0.18 


0.60-0.72 


0.652 


96.3 


2.57 


± 


0.08 


± 


0.22 
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112.0 


0.75 




0.04 




0.09 


0.72-0.90 


0.792 


96.0 


0.91 


± 


0.04 


± 


0.10 


0.792 


111.7 


0.21 


± 


0.02 


± 


0.04 


0.90-1.25 


1.018 


95.9 


0.19 


± 


0.02 


± 


0.04 


1.012 


111.8 


0.03 


± 


0.01 


± 


0.02 



19 



Table 3: Double-differential inclusive cross-section d^a/dpdO, [mb/(GeV/c sr)] of the production of 
tt'^'s in p -I- Be — >^ 7r+ -I- X interactions with -1-8.9 GeV/c beam momentum; the first error is statistical, 
the second systematic; pr in GeV/c, polar angle 6 in degrees. 
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11.83 


± 


0.24 


± 


0.41 


0.30-0.36 


0.329 


66.9 


13.61 


± 


0.26 


± 


0.39 


0.328 


81.9 


8.09 


± 


0.20 


± 


0.26 


0.36-0.42 


0.388 


66.8 


9.72 


± 


0.21 


± 


0.30 


0.388 


81.9 


5.53 


± 


0.16 


± 


0.21 


0.42-0.50 


0.457 


66.7 


7.24 


± 


0.16 


± 


0.27 


0.458 


81.9 


4.43 


± 


0.13 


± 


0.20 


0.50-0.60 


0.546 


66.7 


4.69 


± 


0.12 


± 


0.24 


0.543 


81.9 


2.40 


± 


0.08 


± 


0.15 


0.60-0.72 


0.653 


66.3 


2.51 


± 


0.08 


± 


0.17 


0.654 


81.3 


1.06 


± 


0.05 


± 


0.09 


0.72-0.90 


0.792 


66.4 


0.89 


± 


0.03 


± 


0.09 


0.791 


81.3 


0.27 


± 


0.02 


± 


0.03 


0.90-1.25 


1.010 


66.0 


0.13 


± 


0.01 


± 


0.02 


1 .002 


81.4 


0.03 


± 


0.01 


± 


0.01 








90 <: < w 










10.-) < e 


^ 12 








Pt 


(pt) 


(e) 


d^cr/dpda 


(PT> 


(9) 


dV/dpdn 


0.13-0.16 


0.145 


97.4 


14.32 


± 


0.37 


± 


0.79 


0.145 


114.9 


11.20 


± 


0.27 


± 


0.53 


0.16-0.20 


0.179 


97.3 


13.31 


± 


0.29 


± 


0.54 


0.179 


114.0 


9.75 


± 


0.22 


± 


0.39 


0.20-0.24 


0.219 


97.2 


11.85 


± 


0.28 


± 


0.43 


0.219 


113.7 


7.40 


± 


0.20 


± 


0.29 


0.24-0.30 


0.268 


96.9 


7.55 


± 


0.19 


± 


0.26 


0.267 


113.6 


4.60 


± 


0.13 


± 


0.19 


0.30-0.36 


0.328 


96.8 


5.25 


± 


0.16 


± 


0.20 


0.327 


113.6 


2.73 




0.10 




0.14 


0.36-0.42 


0.387 


96.6 


3.63 


± 


0.14 


± 


0.18 


0.387 


113.2 


1.55 


± 


0.07 


± 


0.10 


0.42-0.50 


0.457 


96.4 


2.10 


± 


0.09 


± 


0.14 


0.455 


112.8 


0.70 


it 


0.04 


± 


0.06 


0.50-0.60 


0.541 


96.2 


1.07 


± 


0.06 


± 


0.09 


0.536 


111.5 


0.22 


± 


0.02 


± 


0.03 


0.60-0.72 


0.650 


95.9 


0.35 


± 


0.03 


± 


0.04 


0.655 


111.2 


0.07 


± 


0.02 


± 


0.02 


0.72-0.90 


0.788 


95.5 


0.09 


± 


0.02 


± 


0.02 

















20 



Table 4: Double-differential inclusive cross-section d^cr/dpdfi [mb/(GeV/c sr)] of the production of 
7r~'s in p + Be — >^ 7r~ + X interactions with +8.9 GeV/c beam momentum; the first error is statistical, 
the second systematic; pt in GeV/c, polar angle 9 in degrees. 









20 < e 


< 30 










30 <e 


< 40 








Pt 


(Pt> 




d^a/dpdn 


(Pt> 


m 




0.10-0.13 


0.116 


24.7 


52.88 


± 


0.99 




3.67 


0.116 


34.7 


37.22 


± 


0.80 


± 


2.54 


0.13-0.16 


0.145 


24.8 


68.86 


± 


1.07 


± 


3.72 


0.145 


34.7 


43.73 


± 


0.81 


± 


2.33 


0.16-0.20 


0.180 


24.7 


73.99 


± 


0.92 


± 


3.23 


0.180 


34.7 


51.74 




0.77 


± 


2.28 


0.20-0.24 


0.220 


24.8 


70.71 


± 


0.89 


it 


2.59 


0.220 


34.6 


50.06 


± 


0.74 


± 


1.85 


0.24-0.30 


0.269 


24.8 


63.44 


± 


0.68 


± 


1.95 


0.269 


34.7 


45.71 


± 


0.57 


± 


1.41 


30-0 36 


0.329 


24.8 


50.19 


± 


0.60 


± 


1.37 


0.329 


34.8 


38.12 


± 


0.52 


± 


1.04 


0.36-0.42 


0.389 


24.8 


39.48 


± 


0.54 


± 


1.10 


0.389 


34.7 


30.08 




0.46 


± 


0.83 


0.42-0.50 


0.458 


24.8 


28.54 


± 


0.39 


± 


0.93 


0.457 


34.7 


20.57 


± 


0.33 


± 


0.66 


0.50-0.60 


0.545 


24.8 


17.69 


± 


0.28 


± 


0.77 


0.545 


34.8 


13.29 


± 


0.23 


± 


0.56 


0.60-0.72 


0.654 


24.9 


9.49 


± 


0.18 


± 


0.56 


0.653 


34.8 


7.13 


± 


0.15 


± 


0.41 


0.72-0.90 
















0.792 


35.0 


2.79 


± 


0.08 


± 


0.22 








40 < 


< 50 










50 <e 


< 60 








Pt 


(pt> 


(S) 


d^cr/dpdn 


(pt> 


("> 


d^a/dpdO, 


0.10-0.13 


0.116 


44.9 


26.16 


± 


0.66 


± 


1.84 
















0.13-0.16 


0.145 


44.8 


32.69 


± 


0.70 


± 


1.74 


0.145 


54.9 


27.00 


± 


0.62 


± 


1.51 


0.16-0.20 


0. 1 80 


44.9 


35.72 


± 


0.63 


± 


1.58 


0.180 


54.8 


30.60 


± 


0.58 


± 


1.38 


0.20-0.24 


0.220 


44.7 


35.26 


± 


0.61 


± 


1.32 


0.220 


54.8 


29.13 


± 


0.57 


± 


1.10 


0.24-0.30 


0.269 


44.7 


31.62 


± 


0.47 


± 


0.98 


0.269 


54.7 


24.82 


± 


0.43 


± 


0.77 


30-0 36 


0.329 


44.7 


27.66 


± 


0.45 


± 


0.76 


0.329 


54.8 


18.64 


± 


0.36 


± 


0.52 


0.36-0.42 


0.389 


44.8 


20.94 




0.38 


± 


0.60 


0.389 


54.6 


14.03 


± 


0.30 


± 


0.42 


0.42-0.50 


0.457 


44.8 


14.98 


± 


0.27 


± 


0.50 


0.457 


54.7 


10.00 


± 


0.22 


± 


0.35 


0.50-0.60 


0.546 


44.7 


8.98 


± 


0.18 


± 


0.40 


0.545 


54.6 


6.11 


± 


0.15 


± 


0.28 


0.60-0.72 


0.653 


44.9 


5.38 


± 


0.13 


± 


0.33 


0.652 


54.7 


3.33 


± 


0.10 


± 


0.21 


0.72-0.90 


0.793 


44.8 


1.79 


± 


0.06 


± 


0.15 


0.793 


54.4 


1.35 


± 


0.05 


± 


0.12 


0.90-1.25 
















1.008 


54.3 


0.23 


± 


0.02 


± 


0.03 








60 < fi 


< 75 










75 < 9 


< 90 








Pt 


(pt> 


{9} 


d^cr/dpdfl 


(Pt> 


(0) 


d^a/dpdft 


0.13-0.16 


0.145 


67.2 


21.82 


± 


0.46 


± 


1.18 


0.145 


82.4 


15.38 


± 


0.37 


± 


0.86 


0.16-0.20 


0.180 


67.2 


22.68 


± 


0.39 




0.95 


0.180 


82.5 


16.85 


± 


0.32 


± 


0.69 


0.20-0.24 


0.219 


67.3 


21.05 


± 


0.38 


± 


0.76 


0.219 


82.3 


15.18 


± 


0.32 


± 


0.55 


0.24-0.30 


0.269 


66.9 


15.96 


± 


0.27 


± 


0.47 


0.268 


82.0 


10.95 




0.22 


± 


0.33 


\J.J\f — JVJ 


329 


66.8 


12.48 


± 


0.23 


± 


0.36 


327 


81.9 


7.40 




0.18 


± 


0.24 


0.36-0.42 


0.389 


66.8 


9.45 


± 


0.20 


± 


0.30 


0.388 


82.0 


4.85 


± 


0.14 


± 


0.19 


0.42-0.50 


0.456 


66.6 


6.04 


± 


0.13 


± 


0.24 


0.457 


81.5 


3.58 


± 


0.11 


± 


0.18 


0.50-0.60 


0.547 


66.7 


3.85 


± 


0.10 


± 


0.20 


0.544 


81.5 


1.82 


± 


0.07 


± 


0.12 


0.60-0.72 


0.651 


66.5 


1.73 


± 


0.06 


± 


0.12 


0.649 


81.3 


0.97 


± 


0.05 


± 


0.09 


0.72-0.90 


0.793 


66.6 


0.69 


± 


0.03 


± 


0.07 


0.782 


81.4 


0.25 


± 


0.02 


± 


0.03 


0.90-1.25 


1.021 


66.3 


0.10 


± 


0.01 


± 


0.02 


1.037 


80.4 


0.04 


± 


0.01 


± 


0.01 








90 < 6» 


< 105 








105 < e 


< 125 






Pt 


(pt> 


(d) 


d'-'a/dpdn 


(pt> 


(d) 


d'^a/dpdn 


0.13-0.16 


0.145 


97.6 


15.03 


± 


0.36 


± 


0.83 


0.145 


114.7 


11.38 


± 


0.26 


± 


0.53 


0.16-0.20 


0.179 


97.3 


14.31 




0.30 


± 


0.57 


0.179 


114.1 


9.10 


± 


0.19 


± 


0.32 


0.20-0.24 


0.219 


97.0 


11.02 


± 


0.25 


± 


0.39 


0.218 


113.9 


6.71 




0.17 




0.24 


0.24-0.30 


0.267 


96.8 


7.66 


± 


0.18 


± 


0.25 


0.266 


113.4 


3.92 


± 


0.11 


± 


0.17 


0.30-0.36 


0.329 


97.1 


4.97 


± 


0.15 


± 


0.20 


0.326 


113.7 


2.45 


± 


0.09 


± 


0.14 


0.36-0.42 


0.387 


96.7 


2.91 


± 


0.11 


± 


0.16 


0.387 


113.2 


1.44 


± 


0.07 


± 


0.11 


0.42-0.50 


0.454 


96.6 


1.88 


± 


0.08 


± 


0.13 


0.454 


112.5 


0.70 


± 


0.04 


± 


0.07 


0.50-0.60 


0.546 


96.7 


0.82 


± 


0.05 


± 


0.08 


0.539 


111.3 


0.22 


± 


0.02 


± 


0.03 


0.60-0.72 


0.649 


95.7 


0.33 


± 


0.03 


± 


0.04 


0.636 


110.1 


0.05 


± 


0.01 


± 


0.02 


0.72-0.90 


0.780 


96.8 


0.08 


± 


0.02 


± 


0.02 
















0.90-1.25 
















1.020 


114.1 


0.04 


± 


0.02 


± 


0.02 
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Table 5: Double-differential inclusive cross-section d'^a/dpdQ [mb/(GeV/c sr)] of the production of 
protons in 7r+ -I- Be ^ p -I- X interactions with +8.9 GeV/c beam momentum; the first error is statistical, 
the second systematic; in GeV/c, polar angle 9 in degrees. 









20 < 


< 30 










30 < e 


< 40 








PT 


(pt> 




d'^a/dpdn 


(pt> 


{9) 


d'^a/dpdU 


0.20-0.24 


0.220 


24.8 


34.44 


± 


0.76 


± 


1.62 
















0.24-0.30 


0.269 


25.1 


36.77 


± 


0.64 


± 


1.54 


0.270 


34.9 


32.90 


± 


0.66 


± 


1.30 


0.30-0.36 


0.329 


25.0 


32.47 


± 


0.58 


± 


1.22 


0.330 


34.9 


30.81 


± 


0.56 


± 


1.07 


0.36-0.42 


0.389 


24.9 


30.62 


± 


0.56 


± 


1.05 


0.389 


34.9 


27.98 


± 


0.54 


± 


0.89 


0.42-0.50 


0.458 


25.0 


28.19 


± 


0.47 


± 


0.90 


0.459 


34.9 


25.72 


± 


0.45 


± 


0.78 


0.50-0.60 


0.547 


24.9 


23.47 


± 


0.37 




0.77 


0.548 


34.9 


20.87 


± 


0.36 


± 


0.65 


0.60-0.72 


0.656 


24.9 


19.42 


± 


0.31 


± 


0.74 


0.657 


34.8 


15.84 


± 


0.28 


± 


0.59 


0.72-0.90 
















0.802 


34.8 


10.53 


± 


0.19 


± 


0.54 








40 < e 


< 50 










50 < 6» 


< 60 








PT 


(pt) 


{9) 




(PT> 


(9) 


d^CT/dpdn 


0.30-0.36 


0.330 


45.0 


30.20 


± 


0.55 


± 


1.00 
















0.36-0.42 


0.389 


45.0 


28.76 


± 


0.55 


± 


0.85 


0.389 


55.0 


26.46 


± 


0.51 


± 


0.75 


0.42-0.50 


0.459 


45.0 


23.87 


± 


0.44 


± 


0.69 


0.458 


54.8 


20.94 


± 


0.40 


± 


0.61 


0.50-0.60 


0.548 


44.9 


17.94 


± 


0.34 


± 


0.58 


0.547 


54.9 


15.89 




0.32 


± 


0.55 


0.60-0.72 


0.656 


44.8 


13.30 


± 


0.27 


± 


0.51 


0.657 


54.9 


11.27 


± 


0.26 


± 


0.48 


0.72-0.90 


0.800 


44.9 


8.01 


± 


0.18 


± 


0.42 


0.799 


54.8 


5.75 


± 


0.15 


± 


0.33 


0.90-1.25 


1.031 


44.7 


2.74 


± 


0.07 


± 


0.22 


1.027 


54.8 


1.69 


± 


0.06 


± 


0.15 








60 < e 


< 75 












75 < 6» 


< 90 








PT 


(pt) 


{9) 


d'^a/dpdn 


(PT> 


{9) 


d^a/dpdn 


0.36-0.42 


0.389 


67.4 


22.14 


± 


0.37 


± 


0.72 
















0.42-0.50 


0.458 


67.3 


20.09 


± 


0.32 


± 


0.58 


0.457 


81.9 


14.52 


± 


0.27 


± 


0.53 


0.50-0.60 


0.547 


67.2 


14.32 


± 


0.25 


± 


0.51 


0.545 


81.8 


9.33 


± 


0.20 




0.42 


0.60-0.72 


0.655 


67.2 


8.69 


± 


0.18 


± 


0.44 


0.654 


81.7 


4.55 


± 


0.13 


± 


0.30 


0.72-0.90 


0.796 


66.7 


4.07 


± 


0.10 


± 


0.29 


0.798 


81.5 


2.00 


± 


0.07 


± 


0.18 


0.90-1.25 


1.026 


66.3 


0.96 


± 


0.04 


± 


0.11 


1.018 


81.5 


0.47 


± 


0.03 


± 


0.06 








90 < 6» 


< 105 








105 < 9 


< 125 






PT 


(PT> 


{9) 


d^o-/dpdf7 


(PT> 


(9) 


d^cr/dpdn 


0.36-0.42 
















0.388 


113.7 


6.70 


± 


0.18 


± 


0.24 


0.42-0.50 


0.456 


96.9 


8.60 


± 


0.21 


± 


0.41 


0.457 


113.5 


4.45 


± 


0.13 


± 


0.20 


0.50-0.60 


0.545 


96.7 


5.17 


± 


0.15 


± 


0.31 


0.544 


112.8 


2.26 


± 


0.09 


± 


0.17 


0.60-0.72 


0.650 


96.6 


2.26 


± 


0.10 


± 


0.18 


0.647 


112.7 


0.71 


± 


0.05 




0.08 


0.72-0.90 


0.792 


96.4 


0.76 


± 


0.05 


± 


0.08 


0.793 


112.1 


0.22 


± 


0.02 


± 


0.04 


0.90-1.25 


1.018 


95.8 


0.19 


± 


0.02 


± 


0.03 


1.016 


111.6 


0.03 


± 


0.01 


± 


0.01 
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Table 6: Double-differential inclusive cross-section d^a/dpdO [mb/(GeV/c sr)] of the production of 
TT+'s in tt"^ -I- Be — >^ tt"*" -I- X interactions with -1-8.9 GeV/c beam momentum; the first error is statistical, 
the second systematic; px in GeV/c, polar angle 9 in degrees. 









20 < 


< ,30 










30 < 


< 40 








PT 


(pt) 




d'^cr/dpdn 


(Pt> 


{d) 


d^o-/dpdt2 


0.10-0.13 


0.115 


24.6 


52.73 


± 


1.28 


± 


3.46 


0.116 


34.7 


38.43 


± 


1.06 


± 


2.58 


0.13-0.16 


0.145 


24.6 


68.39 


± 


1.35 


± 


3.63 


0.146 


34.7 


49.02 


± 


1.16 


± 


2.55 


0.16-0.20 


0.180 


24.6 


85.35 


± 


1.26 


± 


3.81 


0.180 


34.7 


54.25 


± 


0.99 


± 


2.36 


0.20-0.24 


0.220 


24.7 


90.83 


± 


1.28 


± 


3.47 


0.220 


34.7 


59.18 


± 


1.03 


± 


2.21 


0.24-0.30 


0.270 


24.6 


86.94 


± 


1.01 


± 


2.86 


0.270 


34.7 


55.39 


± 


0.80 


± 


1.77 


0.30-0.36 


0.329 


24.7 


78.97 


± 


0.96 


± 


2.38 


0.329 


34.6 


49.51 


± 


0.75 


± 


1.41 


0.36-0.42 


0.389 


24.6 


65.31 


± 


0.86 


± 


1.97 


0.389 


34.5 


41.04 


± 


0.68 


± 


1.14 


0.42-0.50 


0.458 


24.6 


50.82 


± 


0.66 


± 


1.77 


0.458 


34.7 


31.00 


± 


0.51 


± 


0.97 


0.50-0.60 


0.546 


24.6 


34.19 


± 


0.47 




1.65 


0.546 


34.6 


20.89 




0.37 


± 


0.89 


0.60-0.72 


0.653 


24.5 


20.15 


± 


0.32 


± 


1.44 


0.655 


34.6 


12.10 




0.25 


± 


0.76 


0.72-0.90 
















0.796 


34.4 


5.44 


± 


0.13 


± 


0.55 








40 < 


< 50 










50 < 61 


< 60 








Pt 


(Pt> 




d^(7/dpdQ 


(Pt> 


{d) 


d^(T/dpdf2 


0.10-0.13 


0.116 


44.8 


27.23 


± 


0.89 


± 


1.89 
















0.13-0.16 


0.145 


45.0 


34.43 


± 


0.95 


± 


1.80 


0.146 


54.8 


26.57 


± 


0.81 


± 


1.47 


0.16-0.20 


0.180 


44.7 


40.25 


± 


0.86 


± 


1.76 


0.180 


54.8 


29.74 


± 


0.74 


± 


1.31 


0.20-0.24 


0.220 


44.7 


41.13 


± 


0.86 


± 


1.55 


0.220 


54.7 


29.17 


± 


0.72 


± 


1.09 


0.24-0.30 


0.269 


44.8 


37.59 


± 


0.66 


± 


1.21 


0.269 


54.7 


27.17 


± 


0.56 


± 


0.86 


0.30-0.36 


0.329 


44.7 


31.47 


± 


0.60 


± 


0.91 


0.329 


54.6 


22.68 


± 


0.51 


± 


0.66 


0.36-0.42 


0.389 


44.6 


26.25 


± 


0.55 


± 


0.75 


0.388 


54.6 


16.42 


± 


0.42 


± 


0.49 


0.42-0.50 


0.458 


44.6 


19.54 


± 


0.41 


± 


0.62 


0.457 


54.7 


13.11 


± 


0.33 


± 


0.44 


0.50-0.60 


0.545 


44.6 


13.17 


± 


0.29 


± 


0.53 


0.547 


54.6 


8.35 


± 


0.23 




0.36 


0.60-0.72 


0.654 


44.4 


7.69 


± 


0.20 


± 


0.44 


0.654 


54.5 


5.26 


± 


0.17 


± 


0.31 


0.72-0.90 


0.793 


44.5 


3.72 


± 


0.11 


± 


0.32 


0.791 


54.4 


2.28 


± 


0.09 


± 


0.19 


0.90-1.25 
















1.012 


54.3 


0.38 


± 


0.02 


± 


0.06 








60 < e 


< 75 










75 < 61 


< 90 








Pt 


(Pt> 


{6) 


d^(T/dpdn 


(Pt) 


(0> 


d^(T/dpdf2 


0.13-0.16 


0.146 


67.4 


19.92 


± 


0.56 


± 


1.14 


0.146 


82.6 


16.19 


± 


0.51 


± 


1.00 


0.16-0.20 


0.180 


67.2 


21.98 


± 


0.50 


± 


0.97 


0.180 


82.2 


17.37 


± 


0.44 


± 


0.81 


0.20-0.24 


0.219 


66.9 


21.74 


± 


0.50 


± 


0.80 


0.219 


82.2 


16.33 


± 


0.43 


± 


0.59 


0.24-0.30 


0.269 


67.2 


18.47 


± 


0.38 


± 


0.57 


0.268 


81.9 


11.54 


± 


0.29 


± 


0.35 


0.30-0.36 


0.329 


66.8 


14.14 


± 


0.33 


± 


0.41 


0.328 


82.1 


8.72 


± 


0.26 


± 


0.27 


0.36-0.42 


0.389 


66.7 


10.89 




0.28 


± 


0.34 


0.389 


81.9 


5.93 




0.20 




0.22 


0.42-0.50 


0.457 


66.7 


7.65 


± 


0.20 


± 


0.30 


0.458 


81.8 


4.88 




0.17 


± 


0.23 


0.50-0.60 


0.545 


66.6 


5.50 


± 


0.16 


± 


0.28 


0.546 


81.8 


2.80 


± 


0.11 




0.17 


0.60-0.72 


0.654 


66.3 


2.95 


± 


0.11 


± 


0.20 


0.653 


81.7 


1.30 


± 


0.07 


± 


0.11 


0.72-0.90 


0.794 


66.2 


1.03 


± 


0.05 


± 


0.10 


0.785 


81.0 


0.43 


± 


0.03 


± 


0.05 


0.90-1.25 


1.021 


65.7 


0.16 


± 


0.02 


± 


0.03 


1 .007 


81.2 


0.05 


± 


0.01 


± 


0.01 








90 < e 


< iO 


> 








iO.5 < 


< 12 








Pt 


(Pt> 


{6) 


d^CT/dpda 


(Pt> 


{9) 


d^CT/dpda 


0.13-0.16 


0.145 


91. \ 


13.79 


± 


0.46 


± 


0.85 


0.145 


114.8 


10.53 


± 


0.33 


± 


0.66 


0.16-0.20 


0.179 


97.3 


12.92 


± 


0.36 


± 


0.60 


0.179 


114.5 


9.68 


± 


0.27 


± 


0.48 


0.20-0.24 


0.219 


97.2 


11.60 


± 


0.35 


± 


0.42 


0.219 


114.0 


7.56 


± 


0.25 


± 


0.27 


0.24-0.30 


0.267 


97.0 


8.27 


± 


0.25 


± 


0.27 


0.267 


113.5 


5.03 


± 


0.17 


± 


0.19 


0.30-0.36 


0.328 


96.9 


5.50 


± 


0.20 


± 


0.21 


0.327 


113.7 


2.86 


± 


0.13 


± 


0.15 


0.36-0.42 


0.389 


97.1 


4.23 


± 


0.18 


± 


0.21 


0.387 


113.6 


1.71 


± 


0.10 


± 


0.12 


0.42-0.50 


0.457 


96.7 


2.37 


± 


0.11 




0.15 


0.457 


112.4 


0.94 


± 


0.06 


± 


0.08 


0.50-0.60 


0.539 


96.1 


1.34 


± 


0.08 


± 


0.12 


0.541 


111.9 


0.31 


± 


0.03 


± 


0.04 


0.60-0.72 


0.652 


96.3 


0.45 


± 


0.04 


± 


0.05 


0.650 


111.9 


0.12 


± 


0.02 


± 


0.02 


0.72-0.90 


0.794 


95.5 


0.10 


± 


0.02 


± 


0.02 
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Table 7: Double-differential inclusive cross-section d^a/dpdO, [mb/(GeV/c sr)] of the production of 
7r""'s in tt"*" -I- Be — >^ 7r~ -I- X interactions with -1-8.9 GeV/c beam momentum; the first error is statistical, 
the second systematic; pr in GeV/c, polar angle 6 in degrees. 









20 < 61 


< 30 






30 < 6» < 40 


PT 


(pt> 






(pt) 




d^cr/dpdQ 


0.1CM).13 


0.116 


24.8 


46.25 


± 


1.15 


± 


3.10 


0.116 


34.8 


33.42 


± 


0.94 


± 


2.33 


0.13-0.16 


0.146 


24.6 


61.24 


± 


1.26 


± 


3.17 


0.145 


34.6 


36.42 


± 


0.92 


± 


1.93 


0.16-^.20 


0.180 


24.7 


65.33 


± 


1.07 


± 


2.89 


0.180 


34.6 


43.91 


± 


0.88 


± 


1.95 


0.20-0.24 


0.220 


24.6 


70.02 


± 


1.10 


± 


2.61 


0.220 


34.6 


44.85 


± 


0.87 


± 


1.69 


0.24-0.30 


0.269 


24.7 


63.57 


± 


0.85 


± 


2.03 


0.270 


34.7 


42.49 


± 


0.69 


± 


1.35 


0.30-0.36 


0.329 


24.7 


51.65 


± 


0.76 


± 


1.45 


0.329 


34.7 


35.06 


± 


0.62 


± 


0.99 


0.36-0.42 


0.388 


24.7 


42.11 


± 


0.69 


± 


1.20 


0.389 


34.6 


28.74 


± 


0.56 


± 


0.81 


0.42-0.50 


0.457 


24.7 


31.10 


± 


0.51 


± 


1.03 


0.458 


34.7 


19.88 


± 


0.40 


± 


0.65 


0.50-0.60 


0.546 


24.7 


19.63 


± 


0.36 


± 


0.86 


0.546 


34.8 


12.89 


± 


0.29 


± 


0.55 


0.60-0.72 


0.652 


24.6 


10.21 


± 


0.23 


± 


0.61 


0.652 


34.8 


6.87 


± 


0.19 


± 


0.41 


0.72-0.90 
















0.793 


34.7 


2.89 


± 


0.10 


± 


0.24 








40 < 61 


< 50 










50 < 6» 


< 60 








Pt 


(Pt> 


(9) 


d'-'a/dpdn 


(pt) 


(9) 


d^a/dpdn 


0.10-0.13 


0.116 


44.7 


22.32 


± 


0.76 


± 


1.63 
















0.13-0.16 


0.145 


44.8 


26.81 


± 


0.79 


± 


1.48 


0.145 


54.9 


20.72 


± 


0.68 


± 


1.34 


0.16-0.20 


0.180 


44.7 


30.04 


± 


0.72 


± 


1.35 


0.181 


54.8 


23.25 


± 


0.63 


± 


1.07 


0.20-0.24 


0.220 


44.7 


29.53 


± 


0.70 


± 


1.14 


0.220 


54.7 


24.07 


± 


0.64 


± 


0.92 


0.24-0.30 


0.270 


44.7 


27.50 


± 


0.55 


± 


0.88 


0.269 


54.8 


20.77 


± 


0.49 


± 


0.66 


0.30-0.36 


0.329 


44.8 


24.30 


± 


0.52 


± 


0.69 


0.329 


54.6 


16.27 


± 


0.41 


± 


0.49 


0.36-0.42 


0.389 


44.8 


18.89 




0.45 


± 


0.56 


0.388 


54.7 


11.66 


± 


0.34 


± 


0.36 


0.42-0.50 


0.458 


44.7 


13.59 




0.33 




0.47 


0.457 


54.7 


9.07 


± 


0.26 


± 


0.37 


0.50-0.60 


0.544 


44.9 


8.19 


± 


0.22 


± 


0.38 


0.546 


54.8 


5.79 


± 


0.18 


± 


0.30 


0.60-0.72 


0.652 


44.7 


4.95 


± 


0.16 


± 


0.32 


0.653 


54.6 


2.94 


± 


0.12 


± 


0.20 


0.72-0.90 


0.797 


44.7 


1.74 


± 


0.07 


± 


0.16 


0.793 


54.6 


1.29 


± 


0.07 


± 


0.12 


0.90-1.25 
















1.020 


54.4 


0.20 


± 


0.02 


± 


0.03 








60 < 61 


< 75 










75 < 6» 


< 90 








Pt 


(pt> 


(6) 


d'-'a/dpdn 


(pt) 




d^a/dpdn 


0.13-0.16 


0.145 


67.2 


16.30 


± 


0.50 


± 


0.98 


0.145 


82.3 


11.29 


± 


0.39 


± 


0.78 


0.16-0.20 


0.180 


67.4 


17.40 


± 


0.43 


± 


0.81 


0.180 


82.3 


12.37 


± 


0.35 


± 


0.65 


0.20-0.24 


0.220 


67.2 


16.22 


± 


0.42 


± 


0.58 


0.220 


82.4 


11.29 


± 


0.34 


± 


0.48 


0.24-0.30 


0.269 


67.0 


13.86 


± 


0.31 


± 


0.42 


0.268 


82.0 


8.66 


± 


0.24 


± 


0.29 


0.30-0.36 


0.329 


66.9 


10.08 


± 


0.26 


± 


0.31 


0.328 


81.8 


5.91 


± 


0.20 


± 


0.24 


0.36-0.42 


0.388 


66.8 


8.08 


± 


0.23 


± 


0.32 


0.388 


82.0 


4.51 


± 


0.17 


± 


0.20 


0.42-0.50 


0.458 


66.7 


5.41 


± 


0.16 


± 


0.25 


0.456 


81.6 


3.29 


± 


0.13 


± 


0.18 


0.50-0.60 


0.546 


66.7 


3.23 


± 


0.11 


± 


0.19 


0.543 


81.3 


1.79 


± 


0.08 


± 


0.13 


0.60-0.72 


0.654 


66.2 


1.60 


± 


0.07 


± 


0.12 


0.652 


81.2 


0.88 


± 


0.06 


± 


0.08 


0.72-0.90 


0.790 


66.4 


0.74 


± 


0.04 


± 


0.08 


0.786 


81.0 


0.30 




0.02 




0.04 


0.90-1.25 


1.017 


66.2 


0.12 


± 


0.01 


± 


0.02 


1 .004 


81.4 


0.04 


± 


0.01 


± 


0.01 








90 <: < w 










10.-) < e 


^ 12 








Pt 


(pt) 


(e) 


d^a/dpdn 


(PT> 


(6) 


dV/dpdn 


0.13-0.16 


0.145 


97.3 


10.63 


± 


0.38 


± 


0.77 


0.145 


114.9 


8.61 


± 


0.28 


± 


0.72 


0.16-0.20 


0.179 


97.1 


10.19 


± 


0.31 


± 


0.63 


0.179 


114.0 


6.45 


± 


0.20 


± 


0.44 


0.20-0.24 


0.219 


97.0 


7.74 


± 


0.26 


± 


0.43 


0.218 


113.6 


5.65 


± 


0.20 


± 


0.28 


0.24-0.30 


0.267 


97.0 


6.15 


± 


0.20 


± 


0.27 


0.267 


113.8 


3.38 


± 


0.13 


± 


0.16 


0.30-0.36 


0.328 


97.0 


4.37 


± 


0.17 


± 


0.19 


0.326 


113.3 


2.20 


± 


0.11 




0.13 


0.36-0.42 


0.388 


96.7 


2.61 


± 


0.13 


± 


0.16 


0.390 


112.9 


1.27 


± 


0.08 


± 


0.10 


0.42-0.50 


0.457 


96.8 


1.65 




0.09 


± 


0.12 


0.454 


113.6 


0.74 


± 


0.05 


± 


0.07 


0.50-0.60 


0.543 


96.7 


0.89 


± 


0.06 


± 


0.09 


0.540 


111.9 


0.24 


± 


0.03 


± 


0.03 


0.60-0.72 


0.649 


95.8 


0.38 


± 


0.04 


± 


0.05 


0.656 


111.2 


0.09 


± 


0.02 


± 


0.02 


0.72-0.90 


0.791 


96.2 


0.07 


± 


0.02 


± 


0.02 
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Table 8: Double-differential inclusive cross-section d^d/dpdri [mb/(GeV/c sr)] of the production of 
protons in 7r~ -I- Be ^ p -I- X interactions with —8.0 GeV/c beam momentum; the first error is statistical, 
the second systematic; pt in GeV/c, polar angle 9 in degrees. 









20 < e 


< 30 












30 <e 


< 40 








Pt 


(pt) 




d'-'a/dpdn 


(pt> 


{9) 


d'-'a/dpdn 


0.20-0.24 


0.221 


25.1 


28.46 


± 


1.02 


± 


1.39 
















0.24-0.30 


0.269 


25.0 


26.71 


± 


0.79 


± 


1.15 


0.271 


35.0 


26.35 


± 


0.77 


± 


1.07 


0.30-0.36 


0.329 


24.9 


27.13 


± 


0.79 


± 


1.06 


0.330 


34.9 


25.06 


± 


0.75 


± 


0.93 


0.36-0.42 


0.389 


25.0 


24.92 


± 


0.76 


± 


0.89 


0.389 


35.0 


23.28 


± 


0.73 


± 


0.79 


0.42-0.50 


0.459 


25.0 


22.28 


± 


0.60 


± 


0.73 


0.459 


35.0 


17.73 


± 


0.55 


± 


0.57 


0.50-0.60 


0.547 


24.9 


19.55 


± 


0.50 


± 


0.65 


0.548 


34.8 


16.55 


± 


0.48 


± 


0.54 


0.60-0.72 


0.658 


24.9 


14.55 


± 


0.40 


± 


0.57 


0.656 


34.8 


12.50 


± 


0.38 


± 


0.49 


0.72-0.90 
















0.801 


34.9 


7.74 


± 


0.24 


± 


0.40 








40 < 6 


< 50 












50 < 


< 60 








PT 


(pt> 


{9) 


d'-'a/dpdn 


(PT> 


{9) 




0.30-0.36 


0.329 


45.0 


24.35 


± 


0.73 


± 


0.82 
















0.36-0.42 


0.389 


44.9 


20.62 


± 


0.67 


± 


0.63 


0.389 


54.9 


20.67 


± 


0.67 


± 


0.65 


0.42-0.50 


0.457 


45.0 


18.61 


± 


0.56 


± 


0.57 


0.458 


55.0 


16.70 


± 


0.52 


± 


0.53 


0.50-0.60 


0.547 


44.8 


14.12 


± 


0.45 


± 


0.47 


0.548 


55.0 


12.55 


± 


0.43 


± 


0.46 


0.60-0.72 


0.656 


44.9 


9.24 


± 


0.33 


± 


0.37 


0.655 


54.9 


8.68 


± 


0.33 


± 


0.38 


0.72-0.90 


0.799 


45.0 


6.23 


± 


0.23 


± 


0.34 


0.800 


55.0 


4.79 


± 


0.20 


± 


0.28 


0.90-1.25 


1.046 


44.9 


1.78 


± 


0.09 


± 


0.15 


1.044 


54.9 


1.37 


± 


0.08 


± 


0.13 








60 < e 


< 75 














< 90 








PT 


(pt> 


(9) 


d'-'a/dpdn 


(PT> 


{9) 


d-'a/dpdO, 


0.36-0.42 


0.389 


67.5 


18.23 


± 


0.50 


± 


0.61 
















0.42-0.50 


0.458 


67.4 


16.13 




0.42 


± 


0.47 


0.458 


81.7 


10.95 


± 


0.34 


± 


0.42 


0.50-0.60 


0.546 


67.4 


11.55 




0.33 


± 


0.43 


0.546 


81.7 


7.54 


± 


0.26 


± 


0.34 


0.60-0.72 


0.654 


67.2 


6.82 


± 


0.24 


± 


0.36 


0.655 


81.7 


3.88 


± 


0.18 


± 


0.27 


0.72-0.90 


0.797 


67.2 


3.50 


± 


0.15 


± 


0.27 


0.796 


81.7 


1.69 


± 


0.10 


± 


0.15 


0.90-1.25 


1.038 


66.7 


0.78 


± 


0.05 


± 


0.10 


1.024 


81.7 


0.44 


± 


0.04 


± 


0.06 








90 < e 


< 105 








105 < 9 


< 125 






PT 


(pt> 


(9) 


d^a/dpdn 


(pt) 


(9) 


d^a/dpdn 


0.36-0.42 
















0.388 


113.6 


4.91 


± 


0.23 


± 


0.19 


0.42-0.50 


0.457 


96.6 


6.78 


± 


0.27 




0.34 


0.456 


113.6 


3.32 


± 


0.16 


± 


0.16 


0.50-0.60 


0.545 


96.7 


4.23 


± 


0.20 


± 


0.26 


0.541 


113.0 


1.88 


± 


0.12 


± 


0.15 


0.60-0.72 


0.657 


96.5 


2.04 


± 


0.14 


± 


0.17 


0.654 


113.2 


0.79 


± 


0.07 




0.10 


0.72-0.90 


0.795 


96.6 


0.79 


± 


0.07 


± 


0.08 


0.792 


111.9 


0.29 


± 


0.04 


± 


0.05 


0.90-1.25 


1.018 


96.1 


0.23 


± 


0.03 


± 


0.03 


1.035 


111.8 


0.05 


± 


0.02 


± 


0.02 
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Table 9: Double-differential inclusive cross-section d^a/dpdO, [mb/(GeV/c sr)] of the production of 
tt'^'s in 7r~ -I- Be — >^ tt"^ -I- X interactions with —8.0 GeV/c beam momentum; the first error is statistical, 
the second systematic; pr in GeV/c, polar angle 6 in degrees. 









20 < 61 


< ,30 








.30 < 6» < 40 


PT 


(pt> 




d^o-/dpdt2 


(pt) 


(9) 


d'^a/dpdn 


0.1CM).13 


0.116 


24.5 


40.03 


± 


1.60 


± 


2.82 


0.116 


34.7 


27.55 


± 


1.32 


± 


1.90 


0.13-0.16 


0.146 


24.6 


51.88 


± 


1.74 


± 


2.87 


0.145 


34.8 


38.40 


± 


1.48 


± 


2.07 


0.16-^.20 


0.181 


24.5 


58.09 


± 


1.48 


± 


2.61 


0.181 


34.8 


40.43 


± 


1.25 


± 


1.79 


0.20-0.24 


0.220 


24.6 


61.57 


± 


1.52 


± 


2.38 


0.220 


34.6 
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Table 10: Double-differential inclusive cross-section d^cr/dpdO [mb/(GeV/c sr)] of the production of 
7r~"'s in 7r~ -I- Be — >^ 7r~ -I- X interactions with —8.0 GeV/c beam momentum; the first error is statistical, 
the second systematic; px in GeV/c, polar angle 9 in degrees. 
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We refrain from presenting the wealth of cross-section data also in the form of plots. We 
limit ourselves to three representative figures that show the inclusive cross- sections of secondary 
protons, TT+'s, and vr^'s, produced by beams of protons, tt+'s, and vr^'s. We chose a polar- 
angle range that permits a good comparison of our results with published results from other 
experiments. 

Figure 13 illustrates our measurement of the inclusive cross-sections d^a/dpdfi of proton 
and TT^ production by +8.9 GeV/c incoming protons, in the polar angle range 20° < 9 < 30°. 
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Fig. 13: Inclusive cross-sections as a function of charge-signed px of proton and vr^ production by 
+8.9 GeV/c incoming protons, off beryllium nuclei, in the polar-angle range 20° < 9 < 30°. 



Figure 14 illustrates our measurement of the inclusive cross-sections d'^a/dpdil of proton 
and TT^ production by +8.9 GeV/c incoming tt+'s, in the polar angle range 20° < 9 < 30°. 

Figure 15 illustrates our measurement of the inclusive cross-sections d'^a/dpdil of proton 
and TT^ production by —8.0 GeV/c incoming tt^'s, in the polar angle range 20° < 9 < 30°. 
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HARP-CDP jt^ + Be ^ (p,3T^:t") + X 



+8.9 GeV/c 
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Fig. 14: Inclusive cross-sections as a function of charge-signed pt of proton and tt"^ production by 
-1-8.9 GeV/c incoming 7r"'"'s, off beryllium nuclei, in the polar-angle range 20° < 9 < 30°. 
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Fig. 15: Inclusive cross-sections as a function of charge-signed px of proton and vr^ production by 
—8.0 GeV/c incoming vr" 's, off beryllium nuclei, in the polar-angle range 20° < 9 < 30°. 
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9 Comparison with other experimental results 



We compare our +8.9 GeV/c Be cross-sections with those from other experiments and with the 
results obtained by the HARP Collaboration from the same data that we analysed. The data 
of the other experiments are from E802 and E910 which were obtained with somewhat higher 
beam momenta. A more direct comparison with these data will be given in a forthcoming paper 
where we present cross-sections on Be for all beam momenta from 3 to 15 GeV/c [21]. 



9.1 Comparison with E802 results 

Experiment E802 [22] at Brookhaven National Laboratory measured secondary vr+'s in the 
polar-angle range 5° < 6^ < 58° from the interactions of 14.6 GeV/c protons with beryllium 
nuclei. 



Figure 16 shows their published Lorentz-invariant cross-section of 7r+ and tt^ production 
by 14.6 GeV/c protons, in the rapidity range 1.2 < y < 1.4, as a function of mx — rn^, where 
TTT-T denotes the transverse mass. Their data are compared with our results expressed in the same 
units as used by E802. 
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Fig. 16: Comparison of our cross-sections (black circles) of vr^'^ production by +8.9 GeV/c protons 
off beryllium nuclei, with the cross-sections published by the E802 Collaboration for the proton beam 
momentum of 14.6 GeV/c (open circles); the errors are statistical only. 



We note that both experiments agree in suggesting an exponential decrease of the invariant 
cross-section with increasing mx — itIt,, over two orders of magnitude. Unlike the vr^ cross- 
sections, the 7r+ cross-sections at +8.9 and +14.6 GeV/c exhibit nearly the same slope. In the 
comparison of absolute cross-sections, the E802 normalization uncertainty of (10-15)% is to 
be taken into account on top of the beam energy difference. 
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9.2 Comparison with E910 results 

Experiment E910 [23] at Brookhaven National Laboratory measured secondary charged pions 
in the momentum range 0.1-6 GeV/c from the interactions of 12.3 and 17.5 GeV/c protons 
with beryllium nuclei. This experiment used a TPC for the measurement of secondaries, with 
a comfortably large track length of ~1.5 m. With a magnetic field strength of 0.5 T, this large 
track length renders charge identification and proton-pion separation by dE/dx beyond doubt. 
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Fig. 17: Comparison of our cross-sections (black circles) of ir^ production by +8.9 GeV/c protons off 
beryllium nuclei with the cross-sections published by the E910 Collaboration for proton beam momenta 
of 12.3 (open circles) and 17.5 (open triangles) GeV/c; the errors are statistical only. 



Also here, the E910 data are shown as published, and our data are expressed in the same 
units as used by E910. Although the E910 measurements were made with proton beam mo- 
menta of 12.3 and 17.5 GeV/c, respectively, we note the similar 7r+/7r~ ratio between the cross- 
sections from E910 and our cross-sections from a proton beam momentum of 8.9 GeV/c, shown 



in Fig. 17 In the comparison of absolute cross-sections, the E910 normalization uncertainty of 



<5% is to be taken into account on top of the beam energy differences. 
9.3 Comparison with results from the HARP Collaboration 



Figure 18 shows the comparison of our cross-sections of pion production by +8.9 GeV/c pro- 
tons off beryllium nuclei with the results published by the HARP Collaboration [24], in the 
polar-angle range 0.35 < 9 < 0.55 rad. The latter cross-sections are plotted as published, while 
we expressed our cross-sections in the units used by the HARP Collaboration. 

There is a severe discrepancy between our cross-sections and those reported by the HARP 
Collaboration. We note the difference especially of the 7r+ cross-section, and the difference 
in the momentum range. The discrepancy is even more serious as the same data set has been 
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Fig. 18: Comparison of our cross-sections (black circles) of vr^ production by +8.9 GeV/c protons off 
beryllium nuclei with the cross-sections published by the HARP Collaboration (open circles). 



analysed by both groups. For a discussion of the reasons for this discrepancy we refer to the 
Appendix of this paper. 

10 Summary 

From the analysis of data from the HARP large-angle spectrometer (polar angle 9 in the range 
20° < 9 < 125°), double-differential cross-sections d^(j/dpdf2 of the production of secondary 
protons, TT+'s, and tt^'s, have been obtained. The incoming beam particles were +8.9 GeV/c 
protons and pions, and —8.0 GeV/c pions, impinging on a 5% Aabs thick stationary beryllium 
target. The high statistics of the +8.9 GeV/c data allowed us to determine cross-sections of 
K+ and deuteron production, albeit with lower precision. Our cross-sections for 7r+ and tt 
production agree with results from other experiments but disagree with the results of the HARP 
Collaboration that were obtained from the same raw data. 

Acknowledgements 

We are greatly indebted to many technical collaborators whose diligent and hard work made 
the HARP detector a well-functioning instrument. We thank all HARP colleagues who devoted 
time and effort to the design and construction of the detector, to data taking, and to setting up 
the computing and software infrastructure. We express our sincere gratitude to HARP's funding 
agencies for their support. 



32 



Appendix 

The situation that two groups of authors, in this case the 'HARP Collaboration' (referred to as 
'OH' for 'Official HARP') and us, the HARP-CDP group, publish inconsistent results from the 
same raw data, is unusual and unsatisfactory. Naturally, the question arises as to whose results 
can be trusted. 

Since OH have not withdrawn their results despite heavy criticism from independent review 
bodies [9, 10] and from us, we cannot avoid addressing this question. 

The central problem in OH's data analysis is their lack of understanding TPC track distor- 
tions which leads to: 

- a bias of A(1/pt) — 0.3 (GeV/c)~^ in their reconstruction of TPC tracks; in other words, 
their relative bias increases linearly with and attains some 30% at = 1 GeV/c; 
the bias is such that for particles with positive charge px is decreased, while for particles 
with negative charge px is increased; 

- a resolution of (t(1/px) — 0.6 (GeV/c)~^ which is considerably worse than cr(l/px) — 
0.30 (GeV/c)-^ claimed by OH; and 

- a bad overall RPC time-of-flight resolution of 305 ps and an apparent advance of the 
timing signal of protons with respect to that of pions by ~500 ps ('500 ps effect'). 

These three problems, together with a number of additional mistakes [25, 26], have the 
following fatal consequences for cross- sections of secondary hadron production: 

- cross-section spectra of secondary hadrons are distorted especially in regions where cross- 
sections vary strongly with momentum; and 

- protons and pions are partly confused. 

Discussions of the flaws in the OH analysis have been pubUshed in Refs. [3] and [4]. We 
summarize here the main arguments. 

OH's argument that 'dynamic' TPC track distortions can be neglected during the first third 
of the 400 ms long accelerator spill [8], reads as "...owing to their limited mobility the first 
[argon] ions created in the amplification region need about 25 ms to reach the drift region and 
subsequently the steady flow of ions into this region only starts approximately 100 ms after 
the start of the spill, with a gradual transition between these two regimes..." This argument 
is wrong. With an electric field strength of ~1.7 kV/cm the argon ions need less than 1 ms 
for the relevant distance of 1 1 mm. Therefore, dynamic TPC track distortions increase right 
away approximately linearly with time in the spill. (Dynamic track distortions in the HARP 
TPC attain in the r -0 coordinate ~10 mm at the end of the spill, one order of magnitude larger 
than the typical r-0 resolution; the correct algorithms to cope with TPC track distortions are 
described in Ref. [1] and in ample detail in Refs. [27-30].) 

OH's claims in Ref. [8] that "...The constrained fit [which uses the beam point in addition to 
the TPC clusters] is unbiased with respect to the unconstrained fit [which uses the TPC clusters 
only]..." and "...The weight of the vertex constraint compensates very well the distortions..." 
contradict logic. It is impossible that a circle fit of TPC clusters that shifted away from their 
nominal r -0 positions, together with the (undistorted) beam point, yields the same px as a fit of 
the distorted TPC clusters alone. It is equally impossible that a circle fit of the distorted TPC 
clusters together with the undistorted beam point gives an unbiased px estimate. 

OH's justification of their claim of the absence of a bias in their circle fit stems from Fig. 4 
in Ref. [8]: "...In Fig. 4 it is shown that the vertex constraint does not introduce biases for those 
particle trajectories and that the simulation provides an excellent description of the behaviour 
of the resolution function..." . Their Fig. 4 claims to prove the equality of momentum from OH's 
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'constrained' and 'unconstrained fits'. Tliat tliis claim is wrong is evident from the unphysical 
non-Gaussian shape of the shown distribution. Its cause is a mistake in their calculation of the 
r-(f) error of TPC clusters: their is multiplied by a factor cos^ 20 which assigns clusters an 
unphysically large weight depending on how close they are to the singular values = 45°, 
135°, 225° and 315° in the azimuthal angle. (The mathematical intricacies of this mistake are 
explained in Ref. [25].) 

OH never presented evidence that their px resolution during the accelerator spill is indeed 
(t(1/pt) — 0.30 (GeV/c)~^, and that after TPC track distortion corrections their r-0 residuals 
with respect to an unbiased external coordinated system are compatible with zero across the 
whole active TPC volume. 

In their most recent physics publication [24], OH claim "...Corrections that allow use of the 
full statistics to be made, correcting for such [dynamic] distortions, have been developed.. .and 
are fully applied in this analysis. The obtained results are fully compatible within the statistical 
errors and differential systematic uncertainties with those previously published..." . This claimed 
agreement between data from the first third of the spill without distortion correction, with data 
from the full spill with distortion correction, permits the conclusion that OH's full-spill analysis 
is beset by the same flaws as their earlier analysis of data from the first third of the spill. 

Since OH have a biased track momentum, they observe that the RPC timing signal of protons 
is advanced with respect to the RPC timing signal of (relativistic) pions. This '500 ps effect' 
observation led them to conclude in Ref. [34] "...While this is in itself an interesting effect ...it 
prevents the use of the RPCs as a method to verify the reconstructed momentum scale of heavily 
ionizing particles [protons]." As a consequence, they made no use of the powerful particle 
identification capabiUty from RPC time of flight. The exclusive use of dE/dx from the TPC in 
conjunction with a biased track momentum leads to the partial confusion of protons and pions 
in OH's analysis. 

OH's interpretation of the '500 ps effect' is characterized by statements like "...One possible 
explanation is the fluctuation in arrival time of the first cluster of the primary ionization. This 
fluctuation is smaller for heavily ionizing particles [protons]..." in Ref. [7], or "...An order of 
magnitude estimate of the effect given the propagation velocity of electrons in the gas and the 
chamber gap leads to an order of magnitude of a few 100 ps..." in Ref. [34]. This understanding 
of signal generation is wrong. The anode signal is generated by induction. Hence the (fast) 
propagation of electromagnetic waves across the gas gap is relevant and not the arrival at the 
anode of — in comparison — slowly moving electrons. (The correct mechanism of RPC-signal 
generation is described in Ref. [2].) 

We hold that results and conclusions published by OH cannot be trusted. This refers explic- 
itly to their four physics papers [24,31-33], four technical papers [5-8], two Rebuttals [35,36], 
and one Comment [37], published to date. 
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